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Review  of  Hydraulic  Roughness 
Scales  in  the  Fully  Rough  Regime 

A  review  of  predictive  methods  used  to  determine  the  frictional  drag  on  a  rough  surface 
is  presented.  These  methods  utilize  a  wide  range  of  roughness  scales,  including  roughness 
height,  pitch,  density,  and  shape  parameters.  Most  of  these  scales  were  developed  for 
regular  roughness,  limiting  their  applicability  to  predict  the  drag  for  many  engineering 
flows.  A  new  correlation  is  proposed  to  estimate  the  frictional  drag  for  a  surface  covered 
with  three-dimensional,  irregular  roughness  in  the  fldly  rough  regime.  The  correlation 
relies  solely  on  a  measurement  of  the  surface  roughness  profile  and  builds  on  previous 
work  utilizing  moments  of  the  surface  statistics.  A  relationship  is  given  for  the  equivalent 
sandgrain  roughness  height  as  a  function  of  the  root-mean-square  roughness  height  and 
the  skewness  of  the  roughness  probability  density  function.  Boundary  layer  similarity 
scaling  then  allows  the  overall  frictional  drag  coefficient  to  be  determined  as  a  function 
of  the  ratio  of  the  equivalent  sandgrain  roughness  height  to  length  of  the  surface. 
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1  Introduction 

The  most  important  unresolved  issue  regarding  surface  rough¬ 
ness  in  fluids  engineering  practice  is  how  frictional  drag  (for  ex¬ 
ternal  flows)  or  pressure  drop  (for  internal  flows)  relates  to  the 
particular  roughness  topography.  In  other  words,  which  roughness 
length  scales  best  typify  a  surface  hydraulically?  Researchers  have 
been  working  on  this  problem  since  the  early  experimental  work 
of  Nikuradse  [1]  and  Colebrook  [2].  However,  most  previously 
proposed  roughness  scales,  discussed  in  detail  later  in  this  paper, 
are  not  sufficiently  robust  and  are  valid  only  for  regular  roughness 
or  for  a  limited  range  of  roughness  types.  What  is  needed  in  en¬ 
gineering  practice  is  some  means  of  relating  the  rough  surface 
waveform  to  its  frictional  drag.  In  recent  years,  direct  numerical 
simulations  (DNS)  have  further  elucidated  the  understanding  of 
the  underlying  flow  physics  on  rough  walls  (e.g..  Refs.  [3-5]).  To 
date,  however,  DNS  computations  have  been  limited  to  idealized 
roughness  types  at  a  limited  Reynolds  number  range,  making 
them  incapable  of  resolving  issues  regarding  appropriate  scaling 
for  irregular  three-dimensional  roughness.  Therefore,  the  need  for 
accurate  predictive  correlations  in  engineering  practice  remains. 

One  of  the  most  commonly  employed  fluids  engineering  tools 
is  the  Moody  [6]  diagram.  This  useful  diagram  relates  the  pressure 
drop  in  a  pipe  to  the  relative  roughness  (i.e.,  the  ratio  of  the 
roughness  height  to  the  pipe  diameter)  and  the  Reynolds  number. 
Moody  developed  the  diagram  to  be  used  for  naturally  occurring 
roughness  based  on  the  results  of  Colebrook  [2]  for  smooth  and 
rough  pipe  flow.  However,  recent  results  by  Allen  et  al.  [7]  and 
Langelandsvik  et  al.  [8]  show  that  the  Moody  diagram  signifi¬ 
cantly  overestimates  the  pressure  drop  in  the  transitionally  rough 
flow  regime  for  honed  and  commercial  steel  pipes,  respectively. 
This  clearly  indicates  that  the  Colebrook  roughness  function  used 
in  the  formulation  of  the  Moody  diagram  may  not  be  applicable  to 
a  wide  range  of  roughness  of  engineering  interest.  Fortunately, 
from  an  engineering  standpoint,  the  Moody  diagram  likely  gives  a 
conservative  estimate  of  pressure  drop  for  a  given  roughness  in 
the  transitionally  rough  regime. 

A  more  important  issue  regarding  the  Moody  diagram  is  defin- 
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ing  the  appropriate  roughness  length  scale  to  use  as  the  roughness 
height.  The  diagram  was  developed  using  the  equivalent  sandgrain 
roughness  height  ks  [1],  However,  ks  for  a  generic  roughness  can¬ 
not  be  accurately  assigned  a  priori  and  must  be  determined  experi¬ 
mentally.  The  roughness  values  listed  on  the  diagram  for  copper, 
galvanized  steel,  etc.,  were  determined  from  fitting  a  roughness 
height  to  match  a  measured  pressure  drop  from  an  experiment,  or 
in  other  words,  determining  ks.  The  question  then  is  how  the 
height  listed  on  the  diagram  relates  to  a  physical  roughness  length 
scale  measured  from  surface  topography.  This  is  especially  perti¬ 
nent  since  many  of  the  manufactured  surfaces  listed  on  the  Moody 
diagram  do  not  have  a  consistent  roughness.  For  example,  some 
surfaces  are  listed  with  a  roughness  height  that  spans  an  order  of 
magnitude.  Additionally,  the  equivalent  sandgrain  roughness 
height  is  unknown  for  numerous  other  surfaces  of  engineering 
interest.  As  a  result,  one  can  only  expect  to  obtain  accurate  results 
using  the  Moody  diagram  if  ks  is  known  for  the  surface  of  interest 
and  the  flow  is  fully  rough. 

The  naval  architecture  community  includes  the  effects  of  sur¬ 
face  roughness  (e.g.,  paint,  corrosion,  and  fouling)  in  an  allow¬ 
ance  coefficient,  which  is  added  to  the  smooth  surface  friction  and 
residual  (wave  making)  resistance  coefficients  when  determining 
the  overall  drag  of  a  full  scale  ship  [9].  The  International  Towing 
Tank  Committee  [10]  adopted  the  allowance  coefficient  of 
Bowden  and  Davison  [11]  to  be  used  with  the  1978  ITTC  perfor¬ 
mance  prediction  line  for  ship  resistance.  This  coefficient  is  a 
function  of  the  mean  hull  roughness,  the  average  peak-to-trough 
roughness  height  measured  over  50  mm  sampling  lengths  on  the 
hull  surface,  as  compared  with  the  ship  length.  This  equation, 
however,  is  not  an  accurate  hull  roughness  penalty  predictor  since 
it  includes  additional  residual  components  of  resistance  predic¬ 
tion,  including  model  scale  effects  [12],  Townsin  et  al.  [13]  pro¬ 
vided  a  formula  for  predicting  the  roughness  penalty  based  on  the 
mean  hull  roughness  and  the  Reynolds  number.  While  the  inclu¬ 
sion  of  a  Reynolds  number  dependence  allows  for  calculations  in 
the  transitionally  rough  regime,  the  roughness  parameter  is  still 
only  based  on  a  simple  measure  of  the  roughness  height  and  does 
not  account  for  other  roughness  texture  characteristics.  These  con¬ 
siderations,  along  with  a  lack  of  accurate  hull  roughness  measure¬ 
ments,  led  the  ITTC  Specialist  Committee  on  Powering  Perfor¬ 
mance  Prediction  [12]  to  conclude  in  2005  that  the  methods  used 
to  correct  for  hull  roughness  and  fouling  are  of  doubtful  accuracy. 

Computational  fluid  dynamics  (CFD)  has  generally  represented 
rough  surfaces  by  a  smooth  surface  with  modified  boundary  con¬ 
ditions  or  near  wall  equations.  Discrete  roughness  models  include 
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roughness  as  an  additional  drag  term  in  the  near  wall  momentum 
equations.  Using  this  approach,  the  effect  of  roughness  is  confined 
to  the  near  wall  mean  flow,  not  influencing  the  outer  layer  turbu¬ 
lence.  Alternately,  roughness  can  be  accounted  for  in  the  turbu¬ 
lence  by  modifying  the  eddy  viscosity  models  [14],  In  this  ap¬ 
proach,  a  wall  offset  is  included  in  the  mixing  length  model  that 
produces  a  nonzero  eddy  viscosity  at  the  wall.  The  mean  flow  log 
region  extends  to  this  new  origin.  The  amount  of  wall  offset  is  a 
function  of  an  empirically  determined  hydraulic  roughness  length. 
For  the  two  layer  approach,  the  wall  layer  model  is  patched  to  the 
outer  layer  model  by  modifying  the  k  boundary  condition  in  the 
k—  £  model,  and  the  co  boundary  condition  in  the  k-w  model  [15]. 
In  the  fully  rough  regime,  the  proposed  k-s  [16]  and  k-a>  [17] 
models  do  not  require  a  wall  function  since  the  log  layer  extends 
all  the  way  to  the  wall.  The  common  feature  of  all  the  models  is 
an  empirically  determined  term  to  accurately  account  for  the 
roughness.  Predictive  correlations  for  a  wide  range  of  rough  sur¬ 
faces  would  provide  computational  models  the  necessary  rough¬ 
ness  length  scale. 

The  development  of  more  accurate  predictive  correlations  relies 
on  a  robust  database  of  experimental  results  that  have  accurate 
frictional  loss  data  and  detailed  surface  topography  information. 
The  approach  taken  in  the  present  research  is  to  expand  on  previ¬ 
ous  work  that  uses  statistical  moments  of  the  surface  profile,  in¬ 
cluding  a  wider  range  of  three-dimensional  roughness.  Successful 
correlation  of  the  momentum  deficit  due  to  roughness  with  surface 
statistics  will  provide  a  method  of  drag  prediction  based  solely  on 
the  surface  roughness  topography. 


2  Background 

The  velocity  deficit  due  to  roughness  has  been  represented  in  a 
variety  of  ways,  based  on  the  smooth  wall  log-law  profile.  Equa¬ 
tion  (1)  represents  the  effect  of  roughness  by  the  roughness  func¬ 
tion  A  U+  or  the  downward  shift  in  the  log-law  profile 

U+  =  —In  y+  +  B  -  AU+  (1) 

K 


where  k  is  the  von  Karman  constant  and  B  is  the  intercept  for  a 
smooth  wall.  This  can  also  be  recast  into  a  relative  roughness 
form,  Eq.  (2),  where  k  is  a  measure  of  the  roughness  height. 

1  y  1 

U+  =  -In-  +  B  -  \U+  +  —In  k+  (2) 

K  k  K 


Atmospheric  boundary  layer  analyses  generally  represent  rough¬ 
ness  by  a  roughness  length  scale  y0,  which  matches  the  smooth 
wall  log-law  profile,  as  shown  in  Eq.  (3).  In  this  equation,  d  is  the 
zero-plane  displacement,  or  the  distance  below  the  top  of  the 
roughness  where  the  mean  flow  is  zero. 


y-rf\ 
To  / 


(3) 


Equation  (3)  is  only  valid  in  the  fully  rough  regime;  thus  y0  is  an 
alternate  scale  that  is  similar  to  the  equivalent  sandgrain  rough¬ 
ness  height  ks,  forcing  collapse  to  a  single  profile  in  the  fully 
rough  regime. 

Figure  1  shows  typical  rough  wall  boundary  layer  profiles  and 
the  corresponding  roughness  functions.  The  results  presented  are 
from  Flack  et  al.  [18]  for  flows  over  sandgrain  and  mesh  rough¬ 
ness.  The  roughness  heights  for  these  surfaces  are  a  significant 
fraction  of  the  boundary  layer  thickness,  resulting  in  large  rough¬ 
ness  functions.  Similar  roughness  functions  are  obtained  for  the 
mesh  and  sandgrain  surfaces  at  approximately  the  same  unit  Rey¬ 
nolds  number  even  though  the  roughness  heights  differ  substan¬ 
tially.  This  indicates  that  the  roughness  height  alone  is  not  ad¬ 
equate  to  scale  the  momentum  deficit  resulting  from  surface 
roughness. 

If  the  shapes  of  the  mean  profile  in  the  overlap  and  outer  region 
of  the  boundary  layer  are  similar  for  smooth  and  rough  walls,  then 
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Fig.  1  Mean  velocity  profiles  for  mesh  and  sandpaper  sur¬ 
faces  from  Ref.  [18] 


the  roughness  function  can  also  be  expressed  as  the  difference  in 
skin  friction  for  the  smooth  wall  C «■  and  the  rough  wall  C»,  at  the 
same  displacement  thickness  Reynolds  number  (Re^>)  as  proposed 
by  Hama  [19]  and  shown  in 

<4> 

This  allows  roughness  function  results  from  rough  surfaces  ob¬ 
tained  in  the  laboratory  to  be  scaled  up  to  full  scale  using  the  same 
outer  layer  similarity  arguments.  Mean  flow  similarity  can  be 
demonstrated  through  collapse  of  smooth  and  rough  wall  data 
when  plotted  in  a  velocity-defect  form.  Similarity  has  been  experi¬ 
mentally  verified  by  a  number  of  researchers  (e.g..  Refs.  [20-23]) 
including  the  present  authors  for  a  wide  range  of  roughness  types 
[18,24-28],  A  representative  velocity-defect  plot  from  Flack  et  al. 
[18]  for  mesh  and  sandgrain  roughness  is  shown  in  Fig.  2. 

Ideally,  the  roughness  function  could  be  determined  using  a 
correlation  based  on  a  physical  measure  of  the  surface  roughness 
instead  of  having  to  rely  on  laboratory  hydrodynamic  tests.  The 
challenge  is  to  determine  the  proper  scales  obtained  from  surface 
measurements  that  effectively  correlate  with  the  roughness  func¬ 
tion  for  a  wide  range  of  roughness  types.  Figure  3  shows  the 
relationship  between  the  roughness  function  and  the  roughness 
height  for  a  variety  of  surfaces.  For  comparison,  the  uniform 
sandgrain  roughness  results  of  Nikuradse  [1]  are  also  presented. 
Three  flow  regimes  are  observed  in  Fig.  3.  When  the  roughness 
Reynolds  number  k+  is  small,  the  flow  is  hydraulically  smooth 
(i.e.,  AI/+=0).  In  this  case  the  perturbations  generated  by  the 
roughness  elements  are  completely  damped  out  by  the  fluid  vis- 
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Fig.  2  Mean  velocity  profiles  in  velocity-defect  form  from  Ref. 
[18] 
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Fig.  3  Roughness  function  results  for  a  variety  of  rough  surfaces 


cosity,  creating  no  additional  drag.  As  k+  increases,  the  flow  be¬ 
comes  transitionally  rough.  In  the  transitionally  rough  regime,  vis¬ 
cosity  is  no  longer  able  to  damp  out  the  turbulent  eddies  created 
by  the  roughness  elements  and  form  drag  on  the  elements,  as  well 
as  the  viscous  drag,  contributes  to  the  overall  skin  friction.  As  k+ 
increases  further,  the  roughness  function  reaches  a  linear  asymp¬ 
tote.  This  asymptotic  region  at  large  values  of  k+  is  the  fully  rough 
regime.  In  this  regime,  the  skin  friction  coefficient  is  independent 
of  Reynolds  number,  and  form  drag  on  the  roughness  elements  is 
the  dominant  mechanism  responsible  for  the  momentum  deficit. 

In  Fig.  3,  the  peak-to- trough  roughness  height  is  used  as  the 
roughness  scale.  A  common  roughness  scale  in  literature  is  the 
equivalent  sandgrain  roughness  height  ks.  This  is  the  roughness 
height  that  produces  the  same  roughness  function  as  the  uniform 
sandgrain  of  Nikuradse  in  the  fully  rough  regime.  Using  the 
equivalent  sandgrain  roughness  height  as  the  roughness  height  in 
Eq.  (2),  and  the  log-law  intercept  for  a  uniform  sandgrain  surface 


in  relative  roughness  form  «*8.5  (Eq.  (2)),  ks  can  be  determined 
for  a  given  roughness  from  its  roughness  function  AI/+  in  the  fully 
rough  regime  using  the  following  relationship: 

B  -  \U+  +  —In  k  +  =  8.5  (5) 

K 

The  equivalent  sandgrain  roughness  height  is  then  used  as  the 
scale  in  Fig.  4.  As  expected,  all  surfaces  asymptote  to  the  same 
line  in  the  fully  rough  regime.  It  should  be  noted,  however,  that 
ks+  is  unable  to  collapse  the  surfaces  in  the  transitionally  rough 
regime.  Figures  3  and  4  highlight  the  fact  that  a  single  measure  of 
the  roughness  height,  i.e.,  k,  ks,  and  fe,.ms,  is  not  able  to  collapse  the 
roughness  function  results  in  all  regimes  for  a  range  of  surfaces. 

A  number  of  questions  remain  regarding  the  relationship  be¬ 
tween  the  roughness  Reynolds  number  and  the  roughness  function 
for  a  generic  roughness.  The  value  of  k+  (or  ks+)  when  the  surface 


Fig.  4  Roughness  function  results  for  a  variety  of  rough  surfaces  using  ks 
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roughness  ceases  to  be  hydraulically  smooth  has  been  shown  to  be 
a  function  of  the  roughness  type.  Second,  the  shape  of  the  rough¬ 
ness  function  in  the  transitionally  rough  regime  varies  depending 
on  the  roughness  type  and  is  not  known  for  most  surfaces.  Addi¬ 
tionally,  the  value  of  k+  that  defines  the  onset  of  the  fully  rough 
regime  is  unknown  for  most  roughness  types.  The  transitionally 
rough  regime  has  previously  been  defined  as  5  <  k  +  <  70,  based 
on  the  uniform  sandgrain  results  of  Nikuradse  [1],  However,  a 
wide  range  of  values  has  been  reported  in  literature  for  other 
roughness  types.  Ligrani  and  Moffat  [29]  reported  that  the  transi¬ 
tionally  rough  regime  spans  15  <K+<  50  for  a  packed  sphere  bed. 
This  range  is  reported  as  3.5<ks+<30  for  honed  pipe  roughness 
by  Shockling  et  al.  [30]  and  is  given  by  the  present  authors  [26]  as 
2.5<ks+<25  for  a  similar  surface  created  by  surface  scratches. 
Langelandsvik  et  al.  [8]  indicated  that  the  range  of  the  transition- 
ally  rough  regime  is  1.4  <v<  18  for  a  commercial  steel  pipe. 
Lewkowicz  and  Musker  [31]  found  that  the  onset  of  the  fully 
rough  regime  ranged  from  £rms+  =  17  to  40  for  ship-hull  roughness. 
These  disparate  results  for  various  roughness  types  clearly  illus¬ 
trate  the  difficulty  of  identifying  scaling  parameters  that  are  appli¬ 
cable  in  both  the  transitionally  and  fully  rough  regimes.  As  illus¬ 
trated  in  Fig.  4,  collapse  in  the  fully  rough  regime  at  high 
Reynolds  numbers  does  not  ensure  that  the  transitionally  rough 
regime  has  been  properly  captured.  In  fact,  Clauser  [32]  discussed 
the  difficulty  of  finding  a  roughness  scale  in  the  transitionally 
rough  regime  stating  that  some  roughness  types  produce  rough¬ 
ness  functions  with  a  monotonically  changing  slope  while  others 
have  inflection  points. 


3  Previously  Proposed  Roughness  Function 
Correlations 

The  development  of  correlations  for  the  roughness  function  has 
been  an  area  of  active  research  for  many  years.  Correlations  range 
from  simple  models  based  on  roughness  height  and  pitch  to  more 
complicated  relationships  that  include  density  and  shape  param¬ 
eters,  as  detailed  below. 

Bettermann  [33]  proposed  a  functional  relationship  between  the 
rough  wall  log-law  intercept  and  a  roughness  spacing  parameter 
A.,  where  A=pitch/height  of  2D  transverse  bars,  as  shown  in  Eq. 
(6),  for  the  range  of  A  listed. 

/(A)  =  B  -  AU+  +  —In  k+  =  17.35(1.165  log1(l  A  -  1),  1  <A<5 

K 

(6) 

Dvorak  [34]  modified  this  relationship,  using  a  density  parameter, 
A  =  total  surface  area/total  roughness  area,  which  is  equivalent  to 
the  spacing  parameter  of  Bettermann  for  square  bars.  Dvorak 
added  another  relationship,  which  extends  the  range  of  applicabil¬ 
ity  for  more  sparse  roughness  as  follows: 

/(A)  =  17.35(1.165  log10  A  -  1)  1  <  A  <  4.68 

/(A)  =  - 5.95(1.103  log,,,  A- 1)  A  >  4.68  (7) 

These  correlations  were  developed  using  sandgrain  surfaces 
[1,35,36],  mesh  screens  [19],  staggered  rows  of  spheres  [37],  and 
square  bars  [33].  Simpson  [38]  modified  the  parameter  further 
using  A^.— total  surface  area/total  roughness  frontal  area  normal  to 
the  flow,  showing  reasonable  agreement  for  spheres  and  cones 
[39],  staggered  hemispheres  [40],  and  machined  groves  [41]. 

Dirling  [42]  introduced  a  combined  density  and  shape  param¬ 
eter  A.  Roughness  density  is  included  as  the  ratio  of  the  average 
element  spacing  d  to  roughness  height  k,  whereas  shape  is  ac¬ 
counted  for  in  the  frontal  area  of  a  single  roughness  element  Af 
and  the  windward  wetted  surface  area  of  a  single  roughness  ele¬ 
ment  As. 


A  = 


d\  A 


(8) 


Dirling  provided  correlations  for  the  ratio  of  the  equivalent 
sandgrain  roughness  height  ks,  to  roughness  height  k,  for  two- 
dimensional  square  rods,  hemispheres,  spheres,  cones,  right 
angles,  as  well  as  Nikuradse  sand. 


£ 

—  =  0.0164A3,78,  A  <4.93 
k 


-  =  139A~'-90,  A  >  4.93  (9) 

k 

Sigal  and  Danberg  [43]  introduced  the  effect  of  roughness  density 
in  a  different  manner,  including  5,  the  reference  surface  area  be¬ 
fore  adding  roughness,  and  Sf,  the  total  frontal  area  of  the  rough¬ 
ness. 


A  = 


(10) 


The  following  correlations  were  provided  for  two-dimensional 
transverse  roughness  including  bars,  rods,  and  ribs. 


Ic 

—  =  0.00321A  4'925,  1.400  <  A  <4.890 


k. 

—  =  8,  4.890  <  A  <  13.25 
k 


Ic 

—  =  151.711A/1  1379,  13.25  <AS<  100  (11) 

k 

van  Rij  et  al.  [44]  expanded  the  use  of  this  parameter  to  three- 
dimensional  regular  roughness,  using  the  results  of  Schlichting 
[39]  for  staggered  patterns  of  spheres,  spherical  segments,  and 
cones.  The  data  used  incorporated  the  corrections  given  in  Ref. 
[45],  The  correlations  for  three-dimensional  roughness  are  listed 
as  follows: 


—  =  1.583  X  10“5At5'683,  A.  <  7.842 
k  s 

—  =  1.802 A  .003038,  7.842  <  A.  <  28.12 

k 

y  =  255.5AS-1'454,  As  >  28.12  (12) 

van  Rij  et  al.  [44]  also  tested  irregular,  three-dimensional  rough¬ 
ness  and  applied  a  modified  version  of  the  roughness  parameter, 
defined  as 


where  S,  Sf  (as  defined  by  Sigal  and  Danberg  [43]),  and  Ss,  the 
total  windward  wetted  surface  area,  are  calculated  numerically 
based  on  detailed  surface  profiles.  The  equivalent  sandgrain 
roughness  is  determined  from  the  correlations  developed  for 
three-dimensional  regular  roughness,  Eq.  (12),  using  this  modified 
roughness  parameter  and  Aavg,  the  average  roughness  height.  Ex¬ 
perimental  results  for  two  rough  surfaces  were  used  to  test  the 
correlation:  one  completely  covered  by  the  roughness  and  one 
where  rough  regions  of  the  surface  alternate  with  smooth  regions. 
The  predicted  skin  friction,  using  the  calculated  values  of  ks,  com¬ 
pared  well  to  experimentally  results  obtained  from  the  pressure 
drop  in  a  channel. 

Additional  three-dimensional  roughness  predictive  correlations 
were  presented  by  Waigh  and  Kind  [46],  Their  relationships  were 
based  on  the  results  for  a  variety  of  roughness  element  shapes 
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including  cubes,  blocks,  flat  plates,  cylinders,  rods,  cones, 
spheres,  and  hemispheres.  They  expressed  the  effect  of  roughness 
as  a  log-law  shift  C  (Eq.  (14)),  using  a  similar  form  of  the  law  of 
the  wall  as  Eq.  (2). 

I  y  1  1  V 

U+  =  -In-  +  B  -  AU+  +  -In  k+  =  -In-  +  B-C  (14) 

K  k  K  K  k 


The  element  distribution  is  characterized  by  a  density  parameter 
Ak,  which  modifies  the  density  parameter  of  Simpson  [38]  \k, 
with  the  streamwise  aspect  ratio  kl  sm,  where  k  is  the  roughness 
height  and  sm  is  the  streamwise  roughness  length. 


A* 


(15) 


The  following  relationships  were  presented  for  two  density  re¬ 
gimes: 


C  =  10.56  log10 


7.59,  Ak<6 


C  =  -  5.75  log10 


1.38 


+  5.78, 


A,  >6 


(16) 


where  Aw  is  the  roughness  wetted  area,  Ar  is  the  projected  frontal 
area,  and  bm  is  the  spanwise  roughness  length.  Waigh  and  Kind 
discussed  the  fact  that  the  relationships  were  obtained  for  simple 
roughness  element  shapes  in  regular  patterns  and  stressed  the  need 
to  find  ways  to  parametrize  the  geometry  of  natural  roughness. 

Turbine  blade  roughness  that  included  surfaces  with  pitting, 
erosion,  and  deposits  was  studied  by  Bons  [47],  Bons  also  adapted 
the  Sigal  and  Danberg  parameter  for  irregular  roughness  to  calcu¬ 
late  the  equivalent  sandgrain  roughness  for  the  surfaces  as  fol¬ 
lows: 


log(M  =  -  1.31  log(As)  +  2.15  (17) 

The  surface  parameters  were  determined  numerically  from  de¬ 
tailed  surface  profiles,  with  the  roughness  height  k,  taken  as  the 
average  of  the  local  peak-to-trough  roughness  height  for  smaller 
subregions.  The  skin  friction  correlation  using  ks  proposed  by 
Schlichting  [39]  fits  the  experimental  data  best. 

Many  of  the  proposed  correlations  are  valid  for  two-  and  three- 
dimensional  regular  roughnesses,  including  bars,  blocks,  cones, 
hemispheres,  etc.  Determining  the  required  shape  and  density  pa¬ 
rameters  for  an  idealized  roughness  is  a  relatively  straightforward 
procedure.  Extending  this  work  to  naturally  occurring,  irregular 
roughness  is  more  challenging  due  to  potentially  complex  surface 
features.  While  it  is  possible  to  determine  the  required  shape  and 
density  parameters  for  irregular,  three-dimensional  roughness 
[44,47],  complex  numerical  fitting  procedures  are  needed.  A  sim¬ 
pler  method  of  utilizing  the  surface  topography  map  would  be 
useful  for  predicting  frictional  drag  due  to  roughness  in  a  wider 
range  of  engineering  applications. 

Correlations  that  relate  the  statistical  moments  of  the  surface 
profile  with  the  frictional  drag  or  roughness  function  have  been 
offered  by  various  investigators.  Musker  [48]  proposed  an  effec¬ 
tive  roughness  height  to  correlate  the  roughness  function  Af/+ 
=f(k+),  as  shown  in 


(1  +flsp)(l  +bskku) 


(18) 


tive  replicas  of  ship-hull  roughness.  Musker  cautioned  that  the 
cut-off  wavelength  used  to  characterize  the  surface  should  depend 
on  the  size  of  the  roughness  elements  since  the  sizes  of  the  turbu¬ 
lent  eddies  near  the  surface  are  a  function  of  the  longest  surface 
wavelength.  He  proposed  the  use  of  the  Taylor  integral  micro¬ 
scale  [49]  as  the  cut-off  for  the  longest  wavelength  and  two  orders 
of  magnitude  lower  for  the  smallest  wavelengths. 

Medhurst  [50,51]  reported  correlations  for  painted  ship  rough¬ 
ness  based  on  a  parameter  C\h,  which  he  called  the  hydrodynamic 
roughness  number,  having  units  of  roughness  height.  Medhurst 
used  the  following  form  of  the  roughness  function: 


A  U+  =  -In 

K 


C[  hUr 

— - - +  C2 

Bl  v 


(19) 


where  is  found  using  an  alternate  form  of  the  smooth  wall 
log-law  intercept 


B  =  -ln[Bj]  (20) 

K 

Utilizing  the  linearized  form  of  Eq.  (20)  along  with  experimental 
data,  a  regression  analysis  can  be  performed  to  determine  C\h  and 

C2- 

Ble«AU+=(C1h)(^A+BlC2  (21) 


If  the  surface  is  assumed  to  follow  a  Colebrook  form  of  the  rough¬ 
ness  function  (i.e.,  monotonic),  then  C2=1.0  and  A U+  is  only  a 
function  of  Cth.  Medhurst  presented  different  values  of  the  hy¬ 
drodynamic  roughness  number  for  three  types  of  ablative  paints. 
Medhurst  also  cautioned  on  the  use  of  cut-off  wavelengths  and 
other  filtering  techniques,  as  well  as  spatial  resolution  issues  that 
can  bias  results  and  reject  important  surface  information. 

Townsin  and  Dey  [52]  utilized  the  moments  of  the  roughness 
amplitude  energy  spectral  density,  incorporating  a  variable  long 
wavelength  cut-off,  to  determine  the  roughness  function  for 
painted  ship  surfaces.  The  height  of  the  roughness  is  represented 
by  the  spectral  area  m0,  the  slope  by  the  second  moment  m2,  and 
the  range  of  frequencies  composing  the  roughness  profile  by  a 
bandwidth  parameter  a=m0mnl «j22 .  Based  on  the  results  for  26 
ship  paint  surfaces,  the  best  correlation  between  the  roughness 
function  At/+  and  the  roughness  height  k  is  shown  in 

r - 

k  oc  yjam0m 2  =  \  -  (22) 

V  m2 

While  this  technique  shows  promise,  significant  scatter  was 
present  using  this  correlation  in  the  fully  rough  regime  while  the 
correlation  did  not  collapse  the  data  in  the  transitionally  rough 
regime. 

The  correlations  presented  have  demonstrated  some  utility  in 
collapsing  experimental  data  for  a  subset  of  rough  surfaces,  but 
they  yield  a  high  degree  of  variance  when  applied  to  a  larger 
range  of  three-dimensional  roughness  types.  Furthermore,  many 
of  the  correlations  are  difficult  to  apply  to  irregular,  three- 
dimensional  roughness.  The  approach  taken  in  the  present  re¬ 
search  has  been  to  expand  on  the  previous  work  using  statistical 
moments  of  the  surface  profile,  including  a  wider  range  of  three- 
dimensional  roughness.  Successful  correlation  of  the  roughness 
function  with  surface  statistics  will  provide  a  method  of  drag  pre¬ 
diction  based  solely  on  the  surface  roughness. 


4  Full  Scale  Prediction  Methods 


where  krms  is  the  standard  deviation,  sk  is  the  skewness,  ku  is  the 
kurtosis  of  the  surface  elevation  distribution,  and  sp  is  the  average 
slope  of  the  roughness  elements.  The  constants  a  and  b  were 
empirically  chosen  to  collapse  the  data  of  naturally  occurring  sur¬ 
faces.  This  effective  roughness  height  was  reasonably  successful 
in  collapsing  the  roughness  function  for  pipes  covered  with  posi- 


Roughness  function  results  for  rough  surfaces  obtained  in  the 
laboratory  can  be  scaled  up  to  full  scale  using  outer  layer  similar¬ 
ity  arguments.  Thus,  knowing  A U+=f(k+)  allows  for  the  frictional 
drag  of  a  full  scale  rough  surface  to  be  determined.  If  the  equiva¬ 
lent  sandgrain  roughness  ks  is  taken  as  the  roughness  scale,  then 
this  functional  relationship  is  valid  for  all  roughness  in  the  fully 


Journal  of  Fluids  Engineering 


APRIL  2010,  Vol.  132  /  041203-5 


Downloaded  20  Apr  2010  to  131.122.83.98.  Redistribution  subject  to  ASME  license  or  copyright;  see  http://www.asme.org/terms/Terms_Use.cfm 


Fig.  5  Scaled-up  procedure  from  CF  smooth  to  CF  rough 


rough  regime  provided  ks  is  properly  specified.  Schultz  [53]  de¬ 
tailed  the  similarity  methods  used  to  determine  the  overall  fric¬ 
tional  resistance  coefficient  CF  for  rough  wall  boundary  layer  flow 
over  a  flat  plate  of  length  L  if  the  roughness  function  A  U+  is 
known.  The  methodology  incorporates  the  analysis  of  Granville 
[54,55]  and  relies  on  outer  layer  similarity  in  the  mean  flow  for 
smooth  and  rough  walls,  as  demonstrated  in  the  velocity-defect 
profiles  shown  in  Fig.  2.  A  graphical  representation  of  the  scaling 
procedure  is  shown  in  Fig.  5.  Flere,  the  smooth  wall  overall  fric¬ 
tional  drag  coefficient  CF  is  plotted  as  a  function  of  log10(ReL) 
using  the  Karman-Schoenherr  [56]  friction  line  as  follows: 

=  -ln(Re;,  CF)  (23) 

K 

The  rough  surface  overall  frictional  resistance  coefficient  for  a 
known  roughness  function  is  determined  by  displacing  the  smooth 
friction  line  by  a  distance  At/+/c[ln(10)]_1  in  the  positive 
log10(Rei)  direction.  For  a  given  plate  length  L,  a  line  of  constant 
L+=LUjV~l ,  which  satisfies  the  following  relationship,  is  plotted: 


The  intersection  of  this  line  and  the  rough  surface  line  identifies 
CF  for  the  rough  plate  at  a  single  value  of  ReL  for  a  given  A  U+. 

If  this  process  is  repeated  for  a  range  of  Al/+  and  L,  the  overall 
frictional  resistance  coefficient  can  be  mapped  out.  If  the  relation¬ 
ship  A  U+=f(ks+)  is  used  to  account  for  the  roughness,  then  this 
diagram  is  valid  for  all  roughness  in  the  fully  rough  regime  pro¬ 
vided  ks  is  known.  This  is  shown  in  Fig.  6,  as  a  “Moody”  type 
diagram,  where  the  overall  frictional  drag  coefficient  is  presented 
as  function  of  the  ratio  of  ks  to  the  overall  length  of  the  plate  L.  In 
the  fully  rough  regime  the  relationship  can  be  expressed  as  the 
following  formula: 

=  -  2.186  ln(M  +0.495  (25) 

With  an  estimated  10%  uncertainty  in  AI/+  at  95%  confidence,  the 
resulting  uncertainty  in  the  overall  frictional  drag  coefficient  CF  is 
2%  for  ksIL= 0.001,  decreasing  to  less  than  1%  for  ksIL 
=  0.00001. 

On  the  diagram,  it  is  assumed  that  the  onset  of  the  fully  rough 
regime  occurs  at  ks+  =  70,  based  on  the  uniform  sandgrain  results 
of  Nikuradse.  As  noted  previously,  the  onset  of  the  fully  rough 
regime  is  highly  dependent  on  the  specific  roughness  and  likely 
occurs  at  lower  values  of  k*  for  many  types  engineering  rough¬ 
ness.  Additional  measurements  that  span  the  transitionally  rough 
regime  are  needed  to  determine  the  important  scales  for  predicting 
when  a  surface  becomes  fully  rough.  With  the  relationship  given 
in  Eq.  (25),  the  frictional  drag  on  a  planar  surface  can  be  deter¬ 
mined  by  identifying  roughness  scales  that  accurately  predict  ks.  It 
should  also  be  noted  that  the  roughness  function  could  also  be 
incorporated  in  CFD  models  and  used  to  calculate  the  frictional 
drag  of  an  arbitrary  body  covered  with  the  given  roughness. 

5  Development  of  a  New  Roughness  Correlation 

The  development  of  a  roughness  function  correlation  is  re¬ 
stricted  here  to  flows  in  the  fully  rough  regime.  In  this  regime,  a 
larger  number  of  experimental  results  are  available,  and  it  is  pos¬ 
sible  to  collapse  all  roughness  functions  to  a  single  line  if  the 
equivalent  sandgrain  roughness  height  ks  is  used  as  the  roughness 


'°9io  ReL 


Fig.  6  Overall  frictional  drag  coefficient  in  the  fully  rough  regime 
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Table  1  Surfaces  used  in  predictive  correlation  for  ks 


Surface 

Ref. 

k, 

(/im) 

krms 

(/un) 

■s* 

K 

ks  actual 
(H 

ks  predicted 
(H 

%  diff 

Sandpaper — 80  grit 

risi 

546 

67.9 

0.497 

4.49 

529 

522 

1.3 

Sandpaper — 24  grit 

1 18 1 

1291 

167 

0.719 

4.06 

2626 

1954 

25.6 

Sandpaper — 12  grit 

1 18| 

2466 

320 

1.51 

6.17 

6354 

6512 

2.5 

Packed  spheres 

1 28 1 

824 

199 

0.212 

1.90 

876 

1146 

30.8 

Packed  spheres  with  grit 

1281 

738 

158 

0.315 

2.22 

1097 

1018 

7.2 

Honed  (scratch) 

1 26 1 

193 

26.4 

-0.455 

3.63 

71.0 

51.0 

28.2 

Honed  pipe 

L30J 

15.96 

2.5 

0.31 

4.05 

7.5 

16.0 

113 

Commercial  pipe 

rsi 

27.27 

5.0 

-0.19 

2.53 

8.0 

16.5 

106 

Gravel — medium 

I23J 

3079 

605 

0.618 

3.43 

5383 

5167 

4.0 

Gravel — coarse 

[23] 

7350 

1490 

0.0305 

2.46 

6785 

6875 

1.3 

Pyramid — Al 

T27l 

304.8 

72.1 

0.566 

2.40 

510 

589 

15.5 

Pyramid — A2 

127| 

457.2 

108 

0.566 

2.40 

706 

883 

25.1 

Pyramid — A3 

127| 

609.6 

144 

0.566 

2.40 

1301 

1177 

9.5 

Pyramid — B 1 

127) 

304.8 

72.1 

0.566 

2.40 

540 

589 

9.1 

Pyramid — B2 

1 27 1 

457.2 

108 

0.566 

2.40 

577 

883 

53.0 

Pyramid — B3 

[27] 

609.6 

144 

0.566 

2.40 

1012 

1177 

16.3 

scale.  This  scale  understandably  has  its  limitations  since  it  is  not 
intrinsically  related  to  roughness  topography.  The  advantage  of  ks 
is  that  it  provides  a  “common  currency”  among  different  rough¬ 
ness  types  as  pointed  out  by  Bradshaw  [57],  Here  an  investigation 
is  made  into  which  physical  roughness  scales,  if  any,  effectively 
correlate  with  ks.  This  study  was  conducted  using  the  results  from 
a  variety  of  three-dimensional  rough  surfaces,  listed  in  Table  1. 
Numerous  other  researchers  have  measured  the  roughness  func¬ 
tion  in  the  fully  rough  regime  for  other  surfaces.  However,  only 
results  that  also  report  detailed  surface  topographical  measure¬ 
ments  have  been  included. 

A  sample  three-dimensional  topographical  map  for  the  honed 
surface  of  Schultz  and  Flack  [26]  is  shown  in  Fig.  7.  This  surface 
was  profiled  using  a  CyberOptics  laser  diode  point  range  sensor 
laser  profilometer  system.  The  vertical  accuracy  of  the  sensor  is 
1  fim  with  a  laser  spot  diameter  of  10  /rm.  The  data  were  digi¬ 
tized  at  increments  of  25  /rm  in  the  lateral  directions,  and  the 
sampling  area  was  5X5  mm2.  The  other  surfaces  of  the  present 
authors  [18,24,27],  as  well  as  the  surfaces  of  Castro  [23],  were 
profiled  using  the  MicroPhotonics  Nanavea  ST300  white  light 
chromatic  aberration  surface  profilometer.  The  vertical  accuracy 
of  this  system  is  0.3  /rm  with  a  lateral  resolution  of  6  /r m.  Five 
replicate  surface  profiles  were  gathered  on  each  of  these  surfaces 
as  well  as  the  surface  shown  in  Fig.  7.  The  sampling  interval  was 
25  gm,  and  the  sampling  length  was  50  mm.  The  surface  statis¬ 
tics  and  roughness  parameters,  listed  in  Table  1,  were  calculated 


as  the  mean  value  from  the  replicate  profiles.  The  data  were  not 
filtered  to  remove  short  or  long  wavelength  roughness  compo¬ 
nents.  However,  long  wavelength  or  “wavy”  roughness  has  not 
been  included  in  this  study.  Expanding  on  the  work  of  Napoli  et 
al.  [58],  Schultz  and  Flack  [27]  determined  that  roughness  with  an 
effective  slope  ES  less  than  approximately  0.35  does  not  scale  on 
the  roughness  height.  The  relationship  for  the  effective  slope  is 
shown  in 


Filtering,  sampling  interval,  and  sampling  length  can  all  have  a 
significant  effect  on  roughness  statistics  as  pointed  out  by 
Medhurst  [51]  and  Howell  and  Behrends  [59].  This  was  not  a 
primary  focus  in  the  present  study.  The  practical  difficulty  in  ad¬ 
dressing  the  effect  of  these  parameters  is  that  not  all  surfaces  were 
profiled  using  the  same  sampling  interval  and  length.  Specifying 
the  most  appropriate  sampling  interval,  sampling  length,  and  fil¬ 
tering  length  is  quite  complex  and  can  be  expected  to  depend  at 
least  to  some  extent  on  the  inner  and  outer  length  scales  of  the 
flow  in  question.  This  point  underlines  the  need  for  the  adoption 
of  a  consistent  method  of  measuring  hydraulic  surface  roughness, 
as  was  highlighted  by  Howell  and  Behrends  [59]. 

Common  surface  statistical  parameters,  as  well  as  the  wide 
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Fig.  7  Surface  topography  map  and  pdf  of  r,  distance  above  the  mean  roughness  elevation,  for  a  honed  surface 
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Fig.  8  Predicted  ks  versus  actual  ks  using  Eq.  (27) 


range  of  other  roughness  scaling  parameters  reviewed  previously 
in  this  paper,  were  investigated  as  candidate  hydraulic  scales.  This 
was  done  using  a  series  of  statistical  correlations  between  the 
candidate  hydraulic  scales  and  the  measured  equivalent  sandgrain 
roughness  height  ks  for  the  surface.  This  analysis  indicated  that 
the  root-mean-square  roughness  height  scale  (krms)  and  the  skew¬ 
ness  of  the  roughness  surface  elevation  probability  density  func¬ 
tion  (pdf)  ( sk )  had  the  strongest  correlations  with  ks.  Based  on  this 
observation,  it  was  decided  to  develop  a  single  length-scale  cor¬ 
relation  using  both  of  these  parameters.  A  number  of  different 
functional  forms  were  considered.  The  function  that  best  corre¬ 
lates  the  present  data  (Table  1)  is  given  in  Eq.  (27)  and  graphically 
represented  in  Fig.  8. 

K  =f(klms,sk )  »  4.43fems(l  +  a*)1-37  (27) 

In  Fig.  8,  ks  actual  was  determined  using  Eq.  (8)  for  an  experi¬ 
mentally  determined  value  of  A  U+  with  6  =  5.0  and  /c=0.41  in  the 
fully  rough  regime.  The  surfaces  used  in  this  correlation  along 
with  some  of  their  roughness  statistics  are  shown  in  Table  1.  The 
skewness  is  a  quantitative  way  of  describing  whether  the  rough¬ 
ness  has  more  peaks  or  valleys.  A  roughness  with  isolated  large 
peaks  will  have  a  high  positive  skewness.  Surfaces  that  have  be¬ 
come  rough  due  to  deposits  of  roughness  elements  (i.e.,  exhaust 
particulates,  biological  fouling,  etc.)  will  generally  have  positively 
skewed  pdf’s.  A  surface  that  is  pitted  (i.e.,  corrosion,  surface  wear, 
etc.)  will  have  negatively  skewed  pdf’s.  It  should  be  cautioned 
that  only  two  surfaces  used  in  the  present  correlation  had  a  nega¬ 
tive  skewness,  and  both  are  relatively  mild.  The  correlation,  given 
in  Eq.  (27),  has  an  additive  constant  of  1.  This  means  that  the 
present  correlation  would  be  undefined  for  ^<-1.  This  constant 
was  chosen  in  order  to  define  a  predictive  correlation  that  is  near 
the  range  of  the  parameter  space  that  was  investigated.  It  should 
be  noted  that  using  a  larger  additive  constant  would  change  the 
other  constants  in  Eq.  (27)  but  would  not  significantly  reduce  the 
goodness  of  fit.  However,  further  data  for  surfaces  with  negative 
Sj,  are  needed  to  validate  and  refine  this  correlation  for  a  larger 
range  of  the  parameter  space.  Nevertheless,  Fig.  8  shows  the  lin¬ 
ear  correlation  between  the  predicted  value  of  ks  using  the  corre¬ 
lation  and  the  measured  ks.  The  goodness  of  fit  is  excellent  with 
r2  =  0.990. 

Table  1  also  lists  the  difference  between  the  measured  and  pre¬ 
dicted  values  of  the  ks.  The  correlation  works  very  well  for  the 
sandgrain,  sanspray  (gravel),  and  packed  spheres  covered  with 
grit.  This  indicates  that  the  correlation  would  do  an  excellent  job 
of  predicting  the  equivalent  sandgrain  roughness  height  for  large 
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(^>500  /xm)  naturally  occurring  type  roughness.  The  only  ex¬ 
ception  is  the  24  grit  sandpaper  surface.  However,  even  with  a 
26%  difference  in  the  predicted  ks  from  the  measured  value  for  the 
24  grit  sandpaper,  the  predicted  overall  drag  coefficient  CF 
=  0.00773  is  only  6.5%  in  error  from  values  measured  in  tow  tank 
tests  CF=  0.00826  [60]. 

Considering  the  percent  error  from  the  expected  value,  the  cor¬ 
relation  does  not  adequately  predict  the  small  values  of  ks  for  the 
honed  and  commercial  steel  pipes.  This  is  not  surprising  due  to 
the  nature  of  the  fitting  where  the  sum  of  the  square  of  the  residual 
is  being  minimized.  This  will  bias  the  best  fit  toward  the  larger 
roughness  where  the  absolute  value  of  the  difference  is  larger  even 
though  the  relative  difference  is  smaller.  Using  the  results  of 
Shockling  et  al.  [30],  a  113%  difference  in  the  predicted  value  of 
ks  yields  an  error  in  friction  factor  of  approximately  30%.  There¬ 
fore,  the  correlation  presented  may  not  be  adequate  for  all  the 
roughness  considered.  However,  it  appears  that  the  important 
scales  for  better  predictive  correlations  are  the  root-mean-square 
roughness  height  and  the  skewness  of  the  surface  elevation  pdf.  It 
should  be  noted  that  this  is  a  sparse  data  set  to  base  a  predictive 
correlation.  Unfortunately,  experimentally  obtaining  the  rough¬ 
ness  function  for  a  single  roughness  requires  a  large  number  of 
boundary  layer  profile  measurements.  While  pressure  drop  data 
are  easier  to  obtain  in  pipe  flow,  it  is  more  difficult  to  coat / 
manufacture  the  pipe  surface  with  a  specific  roughness.  Recent 
advances  in  computing  the  flow  over  complicated  surfaces,  which 
can  be  easily  parametrically  changed,  may  be  capable  of  provid¬ 
ing  a  larger  data  set  covering  more  of  the  sample  space.  This 
should  allow  for  improvement  in  predictive  capabilities. 

6  Conclusions  and  Recommendations 

An  investigation  has  been  carried  out  to  identify  hydraulically 
relevant  roughness  scales  for  three-dimensional  roughness  in  the 
fully  rough  regime.  A  range  of  common  surface  statistical  param¬ 
eters,  as  well  as  a  host  of  roughness  parameters  from  literature, 
was  considered.  The  results  indicate  that  the  root-mean-square 
roughness  height  (krms)  and  the  skewness  of  the  surface  elevation 
pdf  ( sk )  are  the  most  effective  parameters  in  describing  a  surface 
hydraulically.  A  correlation  based  on  these  parameters  is  offered, 
and  it  shows  promising  agreement  with  the  measured  equivalent 
sandgrain  roughness  height  (ks).  However,  further  data  are  needed 
to  validate  and  refine  this  correlation.  Surfaces  with  negative  sk 
(i.e.,  pitted  or  eroded  surfaces)  would  be  especially  helpful  in  this 
regard.  Collapse  of  the  roughness  function  in  the  fully  rough  re- 
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gime  does  not  ensure  collapse  in  the  transitionally  rough  regime. 
In  fact,  the  nonuniversal  nature  of  the  roughness  function  in  the 
transitionally  rough  regime  makes  a  universal  scaling  for  this  re¬ 
gime  seem  intractable.  Therefore,  the  development  of  a  Moody- 
type  diagram  that  is  applicable  to  a  wide  range  of  surfaces  and  is 
accurate  in  the  transitionally  rough  regime  appears  highly  un¬ 
likely. 
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Nomenclature 


Af  ~ 

As  = 

Ah. 

b  = 

bm  = 

B  = 
C  = 

Cf  = 
CF  = 
C,/i  = 
C2  = 
d  = 
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Fd  = 
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K 
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m0  = 


m  2  = 


mA  — 


N  = 
r  = 
Rea*  = 

ReL  = 

Re„  = 

= 


*p 

S 

Sf 

Sr 


roughness  function  fitting  constant 
frontal  area  of  a  single  roughness  element 
windward  wetted  area  of  a  single  roughness 
element 

wetted  area  of  a  single  roughness  element 
roughness  function  fitting  constant 
spanwise  roughness  length 
smooth  wall  log-law  intercept  **5.0 
log-law  shift  in  a  relative  roughness  form 
skin  friction  coefficient  =  twI \pUe2 
frictional  resistance  coefficient =  FDl\pU2S 
hydrodynamic  roughness  number 
roughness  function  constant 
average  roughness  element  spacing 
effective  roughness  slope 
drag  force 

arbitrary  measure  of  roughness  height 

average  roughness  height 

root-mean-square  roughness 

height  =  \/l  /  AS  J^r2 

equivalent  sandgrain  roughness  height 

maximum  peak-to-trough  height=rmax-rmin 

kurtosis  of  the  roughness  elevation  distribution 

=  (l/A)Ef=1^/[(l/A)Ef=1r?]2 

wavenumber  of  the  surface  roughness 

plate  length 

sampling  length 

zeroth  moment  of  the  roughness  amplitude 
energy  spectra =/q  S, dkx 
second  moment  of  the  roughness  amplitude 
energy  spectra =/gA:2SrtI^_l- 
fourth  moment  of  the  roughness  amplitude  en¬ 
ergy  spectra=/gLxS/ftr 
number  of  samples 

distance  above  the  mean  roughness  elevation 

displacement  thickness  Reynolds 

number  =  UeSt/v 

Reynolds  number  based  on  plate 

length  =  UeLI  v 

momentum  thickness  Reynolds  number  =  Ue6l  v 

skewness  of  the  roughness  elevation 

distribution  =  ( 1  /  N)^=]r]  /  [( 1  /  r2]3/2 

streamwise  roughness  length 

average  roughness  slope 

wetted  surface  area  without  roughness 

total  frontal  area  of  roughness 

energy  spectral  density  of  roughness  elevation 


U  = 
Ue  = 
UT  = 

y  = 

To  = 

a  = 

8  = 
<5*  = 
A  U+  = 

K  — 

X  = 
= 

A  = 
A  k  — 

K  = 

V  = 

e  = 


Superscript 

+  = 

Subscript 

min  = 
max  = 
R  = 
5  = 
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